An atomtronics transistor for quantum gates 
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We present a mechanism for quantum gates where the qubits are encoded in the population 
distribution of two component ultracold atoms trapped in a species-selective triple-well potential. 
The gate operation is a specific application of a new design for an atomtronics transistor where inter- 
species interaction is used to control transport, and can be realized with either individual atoms or 
aggregates like Bose-Einstein condensates (BEC). We demonstrate the operational principle with a 
static external potential, and show feasible implementation with a smooth dynamical potential. 

PACS numbers: 03.67.Lx, 67.85.-d, 05.60.Gg, 03.75.Mn 



As research in ultracold atoms shifts more towards 
practical applications, two of the most promising areas 
are that of quantum computation [l| and atomtronics or 
electronics with trapped ultracold atoms [11-01. Whereas 
classical computation is intimately tied to electronics, 
quantum computation in the context of ultracold atoms 
has so far evolved independently of the relatively new 
field of atomtronics. Taking a cue from classical comput- 
ers, it is likely that analogs of standard electronic compo- 
nents like diodes and transistors can be valuable in quan- 
tum computation as well. In this paper we propose a dif- 
ferent design for an atomtronics transistor which can also 
be used to implement a novel mechanism for a two qubit 
quantum gate. Our proposal has several distinguishing 
features that can be advantageous: (a) qubit encoding 
in the spatial coordinates of particles allows implemen- 
tation with single particle as well as with multi-particle 
entities such as BEC, (b) a system-independent principle 
that offers broad choices for physical realization, (c) op- 
eration does not require manipulation of internal states, 
and (d) easily optimizable for high fidelity and speed. 

A distinguishing feature of our gate mechanism is that 
qubits are directly encoded in and readout from the spa- 
tial distribution of atoms. Spatial mode encoding has 
been primarily used in optical qubits, in the context of 
continuous variable quantum computation [5] , or in dual- 
rail schemes [6]. Certain clever proposals for realizing 
phase gates in double- well potentials 0-0 have also em- 
ployed vibrational modes of trapped atoms. However, 
the gate outcome is contained in the phase of the states 
and readout was shown to require intermediate encod- 
ing on internal states [10]. In contrast, we encode qubits 
in the population distribution of atoms in a triple-well 
potential, and readout simply involves determining the 
presence or absence of atoms in specific wells, possible 
even for single atoms by direct imaging methods [lH. 

Likewise, a simpler operational principle underlies our 
atomtronics transistor, where the atom transport is di- 
rectly controlled by interspecies interaction. Existing de- 
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FIG. 1: Qubit definitions and CNOT quantum gate operation. 



signs are based on manipulating resonant coupling of lat- 
tice sites by adjusting the chemical potential or external 
bias fields (l, y, [HI; or on manipulating atomic internal 
states to transport holes [l3|, or spin (l^ . 

Quantum Gate Operation: We consider a triple- 
well potential (Fig. [T]), with wells labeled left, central, 
rights that can be populated with two repulsively inter- 
acting species A and B. Single qubits are encoded in their 
spatial degrees of freedom: Qubit A is in state |0) or |1) 
if species A is localized in the left well or the right well 
respectively; Qubit B is in state |0) or |1) when species B 
is absent or present in the central well, with 'absent' cor- 
responding to a definite state in an adjacent transversely 
coupled well (see Fig. [3fd)). Our initial objective is to 
create a two-qubit CNOT gate, in a static triple-well, 
such that after a set time T, qubit A is negated if qubit 
B is in |1), but is unchanged if qubit B is in |0). 

Implementation can be understood in terms of the 
three lowest eigenstates (/)o, 0i and 02 for species A in the 
triple- well, with eigenenergies Eq < Ei < E2. As shown 
in Fig. [21 the state 0i has its node in the central well 
while the other two 0o and ^2 have anti- nodes. There- 
fore, when species B is present in the central well, the 
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FIG. 2: The three lowest eigenstates (upper panels) and cor- 
responding eigenenergies (lower panels) of species A in the 
triple- well, with species B absent (left) and present (right). 
The potential is shown with dotted lines. With B absent, 
AE2 = 2 X A^i but with B present, AE^ = 4 x AE^ . The 
energy of the ground state is always the energy reference. 
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repulsive A-B interaction Vab will shift up the energies 
£^0 and £^2, but hardly affect Ei. A class of potentials 
exists where the presence of atom B will raise Eq and 
E2 by the same amount, thus leaving AE2 = E2 — Eq 
unchanged while decreasing AEi = Ei — Eq. 

Species A is prepared in a state {ipAit = 0)) localized in 
one of the two extreme wells. By choosing it to be a Gaus- 
sian with minimized energy for that well, we ensure al- 
most complete projection on the three lowest eigenstates. 
The initial phase relations among the eigenstates, shown 
in Fig. [21 are such that 0o(O) and ^2(0) add up construc- 
tively with (pi (0) in one extreme well and destructively in 
the other. If present, |V^b(0)) is a minimum energy Gaus- 
sian localized in the central well, and we assume a factor- 
ized initial two-particle state \^pAB{^)) = |'^a(O))|'0s(O)) 
(justified at the end of the paper). The process is made 
to work starting from either extreme well by always keep- 
ing the potential bilaterally symmetric about the central 
well minimum. We adjust four degrees of freedom of the 
system: position, width and height of the internal bar- 
riers, and VaB' By simple reparametrization, two are 
used to fix a time scale of operation, and to set the con- 
dition that with species B absent^ the energy differences 
satisfy AE2 = 2 x AEi so that after time T = h/AEi, 
the dynamical phase acquired by the three eigenstates 
are offset by multiples of 27r leading to the revival of the 
initial state in the initially occupied extreme well. The 
remaining two parameters are set to ensure that with 
species B present AE2 = A£^2 remains unaltered while 
AEf = AEi/2 is halved, as seen in Fig. [21 so that now 
AE2 = 4 X AEi . The antisymmetric state ^1 evolves at 
half the rate than without B, hence after the same time 
of evolution T it has an opposite phase, or tt offset, rela- 
tive to the symmetric states. This results in localization 



FIG. 3: Static potential: Evolution of species A in a gate 
cycle with (a) species B absent and (b) species B present, 
(c) The gate frequency 27r/T = AEi/fi varies smoothly with 
Vab. (d) Schematic of transistor, with the central well cou- 
pled to reservoir of species B, and the extreme wells to reser- 
voirs of species A. Transistor Operation: (e) Populations in 
left and right wells after one period t = T, as Vab is varied, 
and (f) the same showing blockade at stronger interaction. 



of the species A in the initially empty extreme well. 

Thus, the presence or absence of the species B leads to 
the revival of species A in the original well or transfer to 
the other extreme well after a set time T, implementing 
the CNOT gate as shown in Fig. [H The evolution of 
the population in the three wells during a gate cycle is 
shown in Fig. [3] (a) and (b) when species A is initially 
localized in the left well. The populations in the wells 
are computed by integrating over the intervals 

(00, —D)^ {—D, D), {D, 00), where ±D are the coordinate 
of the barrier maxima. Figure [3] (c) shows the range of 
gate frequencies (27r/T) for which the gate conditions 
{AE2 = 2 X AEi, AE^ = AE2 and AE^ = AEi/2) can 
be fulfilled and the corresponding interaction strength. 

Transistor Operation: The gate operation described 
above is the specific case, of fixed Vab^ of a general tran- 
sistor mechanism in which species B is used to precisely 
control the flow of species A among the extreme wells. 
A schematic of operation as an atomtronics transistor is 
shown in Fig. [3] (d), where the left and right wells of the 
triple-well setup are coupled to two reservoirs of species 
A, while the central well is coupled transversely via an ad- 
jacent well to a reservoir for species B. Since Vab directly 
affects the dynamic evolution of species A it controls the 
transfer rate from the left well to the right well. Figure [3] 
(c) shows that a range of interaction strengths are avail- 
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able to smoothly adjust AEi which fixes the period for 
a gate cycle T = h/ AEi. Therefore, if we instead fix the 
period of each cycle at T and vary Vab^ the transfer per 
cycle can be controlled; Fig. [3] (e) shows that the frac- 
tion of species A transferred varies smoothly from zero 
to complete transfer. Since Vab depends on the interac- 
tion strength qab and the density IV^bP of species B, for 
large qab small variations in can be used to control 

a large fiux of species A, creating an amplification effect. 
Eventually, as Vab is further increased Fig. [3ff) shows 
that transfer of species A is completely blocked. 

Implementation with trapped atoms: We present 
simulations for ultracold atoms in quasi-one-dimensional 
potentials [15] where only the transverse ground state is 
populated. Species selective potential [Tg] create differ- 
ent directions of relaxed confinement for the two species, 
X for A and y for B, as shown in Fig. [3fd). This keeps 
species B localized (when present) in the central well of 
the triple-well potential along the x-direction, but it is 
free to move in the y-direction, so its presence in the 
central well can be varied. Our analysis can therefore be 
restricted to the dynamics in the x-direction, along which 
we assume harmonic potential for the outer shape of the 
triple- well (excluding reservoirs); it's angular frequency 
00 sets our units: energy hu, length Iq = ^Jlij {mu) and 
time uj~^. The triple- well is completed with two iden- 
tical Gaussian profile barriers symmetrically placed at 
distance, d from the harmonic potential minimum (our 
coordinate origin). For both species the net external po- 
tential along the x-direction has the same form: 

VA{x,t) = \x^^U[e i^+e (1) 

and VB{x^t) = 20 x VA{x^t). In the static model pre- 
sented above Vi{x^t) = Vi{x) is unchanged over time, 
but in the dynamical model we present later, the strength 
U (t) and the position d{t) of the barriers evolve in time. 

The simplest feasible model involves a single atom of 
each species A and B, with effective one-dimensional 
hard-core bosonic interaction strength qab- 

H= \^--g-2^V^{x,,t)j^gABS{xA-XB) (2) 

obtained from the 3D Hamiltonian for hard core bosons 
by integrating out the transverse degrees of freedom [15]. 

The dynamics of species A has little effect on that 
of species B due to its strong confinement, so we may 
factorize the two- atom wavefunction IV^ab) = |'0a)|'0b)- 
By taking the projections (V^bI^IV^ab) and {iPa\H\iPab)^ 
the equations of motion for the Hamiltonian can be re- 
duced to two coupled equations: 

- ^dxxi^A + Va^^A + gABli^B^i^A = idt^A 

-^dxx'ipB + Vb^Pb + 9ab\'iPa\'^iPb = idtipB' (3) 
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FIG. 4: Dynamic Potential: Evolution of species A in a 
gate cycle with (a) species B absent and (b) species B present, 
(c) Populations in left and right wells after one period t = T, 
as Vab is varied, (b) Smooth time evolution of the dynam- 
ical parameters: barrier position d and barrier amplitude U; 
also shown are the energy separations with species B absent 
(A^f , A^f , dotted lines) and present (AEi^AE^). 



The assumption of strong Vb also implies that the in- 
teraction term ^asI'^aPV^b has little effect on the evolu- 
tion of i/jb as we verified. These equations can represent 
two coupled BEC's; self interactions ^a(b) I'^a(b) P can 
be added, but not considered in this paper. We use both 
scenarios for all our results presented here: (i) the two 
particle Hamiltonian in Eq. (|2]) and (ii) the coupled equa- 
tions in Eq. ([3]) and find them almost indistinguishable. 

Dynamic Gate and Transistor: The model pre- 
sented so far assumed a static potential that does not 
vary during operation, and underscores the fact that both 
gate and transistor mechanism operate by evolution of 
the dynamical phases of the states. However, to put this 
in practice, two issues need to be addressed: (i) prepa- 
ration and readout of the quantum bit and (ii) initiation 
and termination of each cycle of operation. Both goals 
can be achieved, provided that the initial and final states 
of species A are the ground state of a single isolated well. 
The population would stay localized in the extreme well 
as long as the barriers are kept high, and the spatial iso- 
lation enables adiabatic loading and readout of the state. 

The gate operation is turned on smoothly by decreas- 
ing barrier height U{t) and separation 2d{t). Being al- 
ready a superposition of the three lowest eigenstates, 
\^|JA{t)) adiabatically follows their evolution, and the evo- 
lution of their relative phases performs the gate operation 
as before. Termination mirrors initiation, since the sec- 
ond half of the process is the reverse of the first half as 
shown in Fig. lU^d). 

Transport algorithm: For simplicity, the path selection 
is based only on the properties of the potential without 
species B; in its presence optimum transfer of species A 
is achieved by a proper choice of the interspecies interac- 
tion. The two parameters U {t) and d{t) are used to ad- 



just AEi and AE2. The barrier width a is kept constant 
here, but could be varied in time as an additional control 
parameter. We start from a potential with high barriers 
and slightly raised bottom of the middle well (with tun- 
neling completely suppressed) we rapidly put the three 
wells into degeneracy, setting AE2 = 2 x AEi = 0.2S?iuj 
by a tiny adjustment in d{t)^ barely visible in Fig. lU^d). 
This step demonstrates the physical isolation of the wells 
achievable at the start and the end of a cycle. 

While maintaining AE2 = 2 x AEi , we slowly increase 
A^i and monitor the phase evolution of the eigenstates 
The speed of variation of the potential is limited by 
the non-adiabatic coupling of the three lowest states to 
the higher order states, defined to be: 
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Due to symmetry of the Hamiltonian, coupling occurs 
only between states of the same parity. The coupling ^02, 
is irrelevant since it simply splits ^0 and E2 further apart 
maintaining AE2 = 2 x A^i. So, we ensure adiabatic 
operation by setting the next higher coupling Ai^ = 0.07, 
a low value. The path selecting algorithm stops when 
the dynamical phase tAEi = tt (half of the revival); the 
second half is obtained by time reversal. 

We ensure that with species B present, exactly the 
same path leads to the desired outcome, of species A lo- 
calizing in the opposite well, by tuning Vab to maximize 
the fidelity. The fidelity is measured by the probability 
of the projection of the final state at t = T on the desired 
CNOT gate outcomes shown in Fig.[TJ and it directly cor- 
relates to the transferred fraction of species A after one 
cycle, plotted in Fig. |4l^c). For both static and dynamic 
potentials, the fidelity was about 98%. 

Discussion and Conclusions: Snapshots [ITJ of the 
dynamics for gate operation shown in Fig.[5]confirm that: 
(i) species B remains localized and is not infiuenced by 
the evolution of species A, (ii) species A is significantly 
delocalized during transit with presence in the central 
well (iii) dynamics with Eq. [21 and Eq. [3] were indis- 
tinguishable on the scale of the plots, suggesting imple- 
mentation with both individual atoms and BECs. With 
BEC, Feshbach resonance can be used to suppress self- 
interaction. With individual atoms optical lattices may 
be preferable and our results would apply with a simple 
reparametrization [l8j. The parameters in our simula- 
tions are U = 13.3huj,a = 2.35/o,<^ = 2.08/o (static) 
and U{0) = 12.9huj,a = 1.61/o,c^(0) = 1.82/o (dynamic) 
and gAB = 0.44 hujlo. Due to proximity of the bar- 
riers in Eq. ([T]), the well depths (barrier-peak to well- 
bottom) never exceed ~ 3huj. On an atom chip [l^] with 

:^ 27r X 5 kHz , the parameters used in our simulation 
correspond to T :^ 1 ms for a gate cycle; similar timescale 
is obtained when we re-parameterize for a lattice. 

After a gate cycle, the reduced density matrix pA{t) = 
TrB|V^AB(t))(V^AB(t)|, gives 1 - Trp^(T) = 4.4 x 10-^ 
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FIG. 5: Snapshots T?] of evolution of population density in 
a gate cycle of period T (time units of cj"^) corresponding 
to Fig. 4(a) and 4(b): upper panels: Density (left axis) of 
species A when species B is absent; the potential is shown as 
blue lines (right axis); lower panels: Joint density of species 
A (horizontal axis) and species B (vertical axis) when species 
B is present. The two-particle Hamiltonian in Eq. ([2]) is used 
here; evolution is identical with the coupled equations Eq. Q. 



indicating a pure state. The von Neumann entropy 
S{T) = -TT[pA{T)ln{pA{T)] = 2.3 x 10-^ with S{t) < 
4.0 X 10~^ during the cycle for the dynamic potential 
and S{t) < 12 X 10"^ for the static. These justify the 
factorization of I^I^ab) and explains the agreement be- 
tween using Eq. [2] and Eq. [3l they also show absence of 
significant entanglement between the species. Gate fi- 
delity can be arbitrarily improved with optimal control 
methods [19] since the mechanism has well-defined ini- 
tial and target states and a highly optimal initial path; 
but we leave out such results as they are system specific. 
This design can implement a universal set of gates, since 
single qubit gates, such as a Hadamard gate, can also be 
implemented by adjusting AEi^ with species B absent, 
for desired partial transfer in a cycle. 
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